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The quantum-kinetic energy of a finite number of trapped fermionic atoms provides a restoring 
force for shear motion due to a distortion of the momentum distribution. In analogy to the twist 
mode of nuclear physics it is proposed that counter-rotating the upper and lower hemisphere of a 
spherical atomic cloud yields a finite-frequency mode closely related to transverse zero-sound waves 
in bulk Fermi liquids. 



The advent of Bose-Einstein condensation of trapped 
atomic 87 Rb in 1995 has initiated large experimental and 
theoretical activities in the field of very dilute, almost 
ideal Bose gases The recent experimental achieve- 

ments of trapping fermionic alkali-atoms J3| raise hope 
that much progress will also be made in the near future 
for Fermi gases. Indeed, Fermi-Dirac degeneracy of a 
mixture of trapped K atoms in two different hypcrfmc 
states has been achieved 0. Recently, also the cool- 
ing of a mixture of bosonic Li and fermionic 6 Li atoms 
into a quantum degenerate regime was accomplished || . 
On the theoretical side, intensive studies have started as 
well. For instance, the possibility that certain gas species 
show attractive interaction (e.g., 6 Li) has initiated stud- 
ies on the possible superfluidity of such systems HQ. 
The close analogy with another finite Fermi system, the 
nucleus, has been pointed out JtJ. It is indeed tempting 
to transpose many typical features of atomic nuclei to 
trapped atomic Fermi gases. Besides the very spectacular 
superfluidity properties there is interest in the spectrum 
of collective excitations, most of them showing features 
proper to Landau's zero-sound modes in bulk Fermi liq- 
uids ||. For finite Fermi systems zero sound translates 
into modes analogous to those of an elastic body pJjo||. 



One of the most remarkable examples is the so-called 
'twist mode' Jll| [l3| in spherical nuclei for which there is 
experimental evidence from backward inelastic electron 
scattering JIJ] . This mode has also been predicted to ex- 
ist in medium to heavy spherical alkali metal clusters as 
the most prominent magnetic multipole excitation [fl5|| . 
In a macroscopic picture the twist mode comprises a co- 
herent oscillation of the particles in the upper half sphere 
versus those in the lower half sphere. For small ampli- 
tudes it corresponds to a purely kinetic excitation with- 
out any spatial distortion of the equilibrium shape. In 
this note we wish to investigate to what extent the twist 
mode may also occur as a collective mode in very dilute, 
atomic Fermi gases in the degenerate limit. 

It is well known that the twist mode is generated by 
the operator Q 

T = e~ iazl * = e aU ^, u = (yz,-xz,0) (1) 

acting on the ground-state wave function of the Fermi 
system. As is evident, T induces a rotation of the par- 
ticles around the body-fixed z-axis with a rotation angle 
proportional to z, i.e., the rotation is clockwise for z > 
and counterclockwise for z < 0. The amplitude of this 
twist is characterized by the angle a. One can verify that 
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the twist corresponds to a magnetic mode of spin-parity 
= 2~ . Although the operator (0) induces no change 
in the spatial distribution, the momenta become locally 
distorted. The subsequent derivation of the mode fre- 
quency will closely follow the original work of Holzwarth 
and Eckart fll2|| . For atomic Fermi gases with N ~ 10 5 
to 10 6 particles the Thomas-Fermi approach is very ap- 
propriate j?J (as is the case of atomic bosons 16 1). Most 
relevant for our purpose is the total kinetic energy of the 
system: 



. , d 3 rd 3 p p 2 

E km {a) = / % ^ —f a {r,p) 



(2) 



where f a is the distorted phase-space distribution in 
Thomas-Fermi approximation 



f a = v8{p 2 F {r 1 p)-p 2 ) 



(3) 



Here v is the degeneracy factor and 9 the unit step 
function. The (quadrupole)deformed local momentum 
is given by 



PF(r,p) = p F (r)N(a) {l 



a^-[y(Y 21 (p)-Y 2 ^(p)) 
+ix(Y 21 (p) + Y 2 _ 1 (p))]} 



(4) 



where 



p F (f) = ^2m[p-V ex {r)-g{v-l)p{r)} (5) 

denotes the local Fermi momentum with chemical poten- 
tial p, and trapping potential V ex (f) and 



N(a) = 1 - a 



2 x 2 + y 2 
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(G) 



In the dilute gas limit it can be assumed that the two- 
body interaction is given by its long-wavelength limit 



v(r — r') = g5(r — r') 



(7) 



with g = Anh a/m, where a is the s-wave scattering 
length. For trapped atomic Fermions two situations can 
occur. For example 40 K can be trapped as a mixture 
of atoms in two different m F states, mp — 9/2 and 
nip = 7/2 0| . In this case s-wave scattering is real- 
ized and the interaction (Q) is active. In contrast, if 
the atoms are all in a single hyperfine state, there can 
be no s-wave scattering. In the latter case p-wave in- 
teractions may become very important |17] . The inclu- 
sion of p-wave interactions as well as the description of 
Fermion-Boson mixtures requires a substantial extension 
of the fluid-dynamical formalism and will be discussed in 
a following publication. In the present work we concen- 
trate on a two component Fermi gas with equal number 
of particles in each magnetic substate. In this case the 



equilibrium Thomas-Fermi equation for the density p(r) 
of atoms in one of the two mp states is given by [|17j|7|l: 



p(r) = 



6ir 2 h 3 



(8) 



which leads to a cubic equation for p which can be solved 
analytically. To second order in a we then obtain 

E km {a) = J d 3 r T (f) [1 + y (x 2 + y 2 )} (9) 

where To is the total kinetic energy density at equilibrium 
3 K 2 p F {r) 



r (r) = 2 



5 2m 6n 2 h 
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(10) 



Since in lowest order, the potential energy contains no 
contribution from the twist mode, we obtain for the 
restoring force (assuming a spherically symmetric trap 
of harmonic oscillator shape with frequency oj) : 



C 



d 2 E km {a) _ 16tt f, 4 



9a 2 



dr r To(r) 



(11) 



We also need to evaluate the mass parameter B of the 
twist motion. As usual in fluid dynamics it is given by 



B = m J d 3 r 2p{r) u 2 



m J d 3 r 2p{r) z 2 {x 2 + y 2 ) 
16tt 
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m I dr r p{r) , 



(12) 




where 2p(r) corresponds to the total density of atoms. 
The twist frequency fix is then obtained as 



(13) 



We have considered two systems. One is 6 Li with a 
very large attractive scattering length of a=-2063ao (ao 
= Bohr radius) [^8). The trapping potential was taken 
to be Tiui = 6.9 nK ||. The other system is 40 K with a 
repulsive scattering length of a=157ao and hw—1.6 nK 
|19| . The results for the twist frequency as a function of 
the particle number in each magnetic substate are given 
in Table 1 . In order to see how Qp depends on the inter- 
action strength (which may be variable due to the tuning 
of Feshbach resonances) we also list Ot for various other 
values of the scattering lengths (differing from the origi- 
nal ones by powers of 10). 

From Table 1 it can be inferred that the influence of 
the interaction on the twist frequency is very moderate. 
It is typically of the order of 10% (except for very large 
particle number). This is consistent with the expecta- 
tion that, for transverse zero sound, s-wave interactions 
give no contribution to the restoring force (l3). The de- 
pendence of the twist frequency on the interaction only 
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enters through the mass parameter B which depends on 
the density (Eq. |l^) Depending on the sign of the in- 
teraction the gas either expands (repulsive) or contracts 
(attractive) relative to the free gas case. Consequently 
the frequency is decreased (increased) with respect to its 
non-interacting value £Ito = to for repulsion (attraction) . 
This feature should be measurable even though the ab- 
solute effect might be small. 

Very interesting possibilities arise when the p-wave in- 
teraction becomes strong jL7|. In this case there will be 
a significant correction to the kinetic energy through the 
effective mass and hence a large influence on the twist fre- 
quency. One may even encounter instabilities, signaled 
by the exponential growth of the twist amplitude. Also, 
as mentioned above, one can then have a direct influ- 
ence of the interaction in a one component Fermi sys- 
tem. Another interesting issue is how the twist mode 
is influenced by eventual superfluidity. It can be pre- 
dicted that the twist mode ceases to exist once pairing 
is strong enough for the system to reach its irrotational 
flow limit Q. There may, however, be intermediate sit- 
uations. To our knowledge these possibilities have not 
been addressed for the nuclear twist. 

The question how to excite the twist mode in the exper- 
iment may not be trivial. One could imagine utilizing the 
well-developed technique of rotating trapped atoms [po| . 
First a very elongated trap potential is created. Subse- 
quently a rotating laser field is wrapped around the long 
axis inducing a rotation of the atomic cloud. If, instead 
of applying the laser field parallel to the long axis, it is 
incident at a certain angle <f> and a mirror is placed paral- 
lel to the long axis which reflects the laser beam in such 
a way that it hits the other hemisphere at an angle — 0, 
then the first hemisphere will rotate in one direction and 
the other hemisphere in the opposite direction (at least 
approximately) since the rotational sense of the laser is 
inverted by the mirror. If the rotation is very gentle and 
stopped at a certain time, the system will continue oscil- 
lating (approximately) in the twist mode. Whether one 
can detect the twist mode by switching off the trap poten- 
tial and subsequently imaging the velocity distribution of 
the expanding atoms remains an open question. 
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TABLE 1: The twist- mode frequencies £It in units of the trap frequency ui for several numbers of 6 Li (g < 0) and 
40 K (g > 0) atoms in each hyperfme state for different scattering lengths (in units of the Bohr radius). The frequencies 
of the harmonic oscillator traps are u> = 2tt x 144 Hz for 6 Li and u> — 2tt x 33.5 Hz for 40 K, which corresponds to level 
spacings of hui = 6.9 nK and 1.6 nK, respectively. 
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